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' Abstract. The stellar line-of-sight velocity distribution (LOSVD) can be strongly asymmetric 
in regions where the light contributions of both disc and bulge in spiral and lenticular galaxies are 
comparable. Existing techniques for the stellar kinematics analysis do not take into account the 
difference of disc and bulge stellar populations. Here we present a novel approach to the analysis 
of stellar kinematics and stellar populations. We use a two-component model of spectra where 
' different stellar population components are convolved with pure Gaussian LOSVDs. For this 
| model we present Monte-Carlo simulations demonstrating degeneracies between the parameters. 
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r~| | 1. Introduction 

£pHj A fiat rotating stellar disc and a slowly rotating spheroidal bulge in spiral and lenticular 
galaxies usually possess very different stellar population properties. Consequently the 
resulting stellar line-of-sight velocity distribution (LOSVD) can be strongly asymmetric 
in regions where the light contributions of both disc and bulge are comparable. At the 
first approximation, this can be accounted by the Gauss-Hermite parametrization of 
the LOSVD (van der Marel & Franx, 1993), however, different absorption features (e.g. 
sensitive to age and metallicity) will have different effective LOSVDs. The first attempt 
to recover parametrically the multi-component dynamics connected to multi-component 
^vq ' stellar populations was done by De Bryuine et al. 2004 where the authors used stellar 
spectra to model different stellar populations. Here we present a more realistic approach 
lO ' based on the full spectral fitting, where multiple dynamical components are represented 
. by multiple simple stellar population (SSP) models. More sophisticated stellar population 
t— I ' models can be used instead of SSPs. 



2. Two-component parametric LOSVD recovery 

As an example, here we present a two-component model having two pure Gaussian 
LOSVD components with different stellar populations characterised by their ages and 
metallicities. An optimal template is represented by the linear combination of two SSPs 
each convolved with its own LOSVD, hence the x 2 value is computed as follows: 

2 _ ^ [F^ P p ■ £ . k 2 ■ SjT^Z,) g C(v 3 , q,)] 2 

where C(v, a) - pure Gaussian LOSVD; Fi and 8F, are observed flux and its uncertainty; 
S(Tj,Zj) is the flux from the j-th synthetic spectrum of SSP with given age Tj and 
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Table 1. Parameter of stellar components. 



Component 


v, km/s 


a, km/s 


Age, Gyr 


Z, dex Weight Kj 


1 (Bulge) 





300 


15 


-0.2 0.7 


2 (Disc) 


100 


100 


5 


0.0 0.3 



metallicity Zj\ P p is multiplicative Legendre polynomials of order p for correcting the 
continuum which determined at each step of minimization loop by solving the linear 
least-square problem; kj is the j-th component weight (normally found by the linear 
minimization). The important point in this study is that we fixed the relative SSP con- 
tributions kj to the values derived from the photometric light profile decomposition. 
This approach was implemented on top of the NBursts full spectral fitting technique 
dChilingarian et al. 2007a| |Chilingarian et al. 2007b[ ). 

Our new approach was applied to the real spectra of luminous early-type galaxy 
NGC 524 (see Katkov et al. (2011)). It was shown that LOSVD of NGC 524 exhibits 
strong asymmetry and a stellar component of galaxy can be described as a bulge, inner 
disc and a conterrotating outer disc. 

3. Monte-Carlo simulation 

In order to explore the parameter space, we constructed and analysed using our tech- 
nique a sample of 8000 realizations of synthetical spectra created as a linear combination 
of two components (see Table [1]) adding the random noise corresponding to the signal- 
noise ratio (SNR) 100. The relative contributions fci,2 were fixed. The degeneracies be- 
tween all the parameters are shown in Fig. Q] as a covariance matrix. Here and after red 
and blue points correspond to the bulge and disc components respectively. We can see 
that the most degenerated pairs of parameters are Ti — T 2 and Z\ — Z 2 ■ 

In order to demonstrate how the parameters are recovered depending on the SNR we 
fitted a sample of 8000 models varying a SNR between 30 and 250. Fig.^top line) shows 
how the output parameter uncertainty decreases when increasing SNR. The dashed lines 
show the input parameters of synthetic spectrum. 

We used the same approach to demonstrate how important is the knowledge of relative 
contributions kj . Fig. [^bottom line) presents 4000 MC realizations where we varied the 
relative SSP contribution of bulge component k\ in range between 0.5 and 0.9 while the 
real value was equal to 0.7. SNR=150 was adopted in these realisations. It is obvious 
that even little offset in kj can lead to important biases in the recovered parameters. 

Authors thank the IAU for the provided financial aid and RFBR grant 10-02-00062 
for covering the remaining expanses. 
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Figure 1. Covariance matrix of the model parameters. Left-top quadrant and right-bottom 
quadrant present a self-correlation of parameters of bulge and disc components, correspondingly. 
Left-bottom and right-top quadrants correspond to mutual correlation of stellar components. 
White crosses show the input values. 



Velocity Dispersion Age Metallicity 




0.5 0.6 0.7 0.8 0.9 0.5 0.6 0.7 0.8 0.9 0.5 0.6 0.7 0.8 0.9 0.5 0.6 0.7 0.8 0.9 
weight 1 weight 1 weight 1 weight 1 



Figure 2. Top line: output parameter dispersion as a dependence on signal-to-noise ratio. 
Bottom line: output parameters with varied relative contribution of components. 



